Very Narrow Resonances in Spherical Proton Emitting Nuclei 
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The separation energy and half-life of some heavy proton emitting nuclei, and the 
single-particle structure of the unbound 11 N, have been evaluated by implementing a 
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tion context. The basic scheme behind the method consists in using the ground-state 



careful numerical treatment to solve Schrodinger equation in a continuum discretiza- 
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. proton emitter in connection with an isolated single-particle resonance. 
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Recently, several unstable proton emitting nuclei with medium and heavy masses fl 



111 ] have been discovered and are attracting much attention in both 
theoretical and experimental nuclear physics. In the light mass limit unbound 11 N is another 



well studied proton emitter nucleus which presents some intriguing phenomena such as the 
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S1/2 intruder level and the u Be mirror states discussed in the recent literature 
The parent nucleus decays by proton emission in a quantum tunneling process. In a first 
approximation we can treat this problem as an unbound proton + core system in which the 
ground-states instabilities are studied from the single-particle resonance point of view. These 
resonances have large half-lives as long as ~ Ins — Is. One of the great difficulties presented 
in these calculations is the determination of the separation energies and half-lives because 
these states have a long lifetime in comparison with the oscillation time inside the potential 
well: the energies are of the order of ~ lMeV and are extremely narrow resonances (T < 
10~ 12 MeV) corresponding to the half-lives ti/2 > Ins (an enormous numerical computation 
time would be necessary to sweep some MeVs of energy by steps of the order of their width 
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size in search of the resonance position ). The spin-parity assignment to the measured levels 
is another arduous task, and because the proton decay has a strong dependence on the 
angular momentum, the spin-parity predictions for some ground-state proton emitters are 
assigned by matching the experimental half-life with the theoretical calculations such as the 
WKB method|5, 81. Several theoretical methods have already been presented 

to calculate the half-lives and spins of these nuclei and we are presenting another method 
based on solving the Schrodinger equation in a continuum discretization context Q|. 
The purpose of this work is to calculate the decay properties, such as energies and half-lives, 
of some proton emitters implementing a careful numerical treatment to this continuum 
projection method. 

The ground-state proton emitters are treated in connection with an isolated long lifetime 
single-particle resonance. Using orthogonal and complementary projectors the scattering 
solution \if) + ) of the single-particle Hamiltonian H at energy E can be written asQ] 

\^ + ) = \u)(u\^ + ) + P\^ + ), (1) 

where \u) is a normalized single-particle state and P — 1 — projects onto the comple- 

mentary sub-space of the single-particle Hilbert space. Projection and formal manipulation 
of the Schrodinger equation (E — H)\tp + ) = then gives 

P|V + > = + G + pp (E)PH\u)(u\4,+), (2) 



«IV> + > = „ ' ':' + ,„ • (3) 
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where Gp p (E) = (E + ir/ — H pp )~ l is the Green's function and |<^ + ) is the scattering solution 
of the projected Hamiltonian H pp at energy E. The \ solution is obtained by solving the 
inhomogeneous equation 

(E-H)\<p + )=a\u), (4) 

with a = — (u\H\lp + ) adjusted to satisfy the orthogonal relation (u\ip + ) = 0. The resonant 
energy, E Q , is defined as the energy for which the integrated internal probability density is 
largest and \u) is chosen as being proportional to the internal part of the resonant wave func- 
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tion \iP + )e=e , truncated at a radius somewhat larger than the potential radius, containing 
the principal behaviour of the resonance inside the well: 

(r\u) = N(l + e^-^r 1 W + >£=£ , (5) 
where N is a normalization factor, R u and a u are the truncation parameters for (r\u). Thus 
the complex single-particle resonance energy (e u — ^f- = (u\H\u) + (u\HGp P (E)H\u)) is 
calculated in the energy E . The imaginary part (T u = 27r|a| 2 ) corresponds to the single- 
particle width and the associated half-life is given as ti/2 = Jr- In 2. The potential utilized in 
the single-particle Hamiltonian contains the central potential in a Woods-Saxon form and 
the spin-orbit term with a Woods-Saxon derivative form, also including the centrifugal and 
Coulomb ( uniformly charged sphere) terms. 

To determine the resonance parameters for the ground-state proton emitters, special at- 
tention must be taken in the numerical implementation of this method because these very 
narrow resonances are difficult to detect and require a more accurate search. The single- 
particle resonance energy e u is obtained imposing that it must reproduce the experimental 
Q value by adjusting the central and spin-orbit potentials. To avoid wasting much compu- 
tational time, and to achieve a better accuracy in the results, we start with an initial energy 
range near the experimental Q value to solve the Schrodinger equation. The energy range 
is swept with a integration energy step AE to localize the energy E where e u and T u are 
calculated. After this, a new energy range is defined near e u . The range size is fixed at 
~ 2 x 10 3 Ai? centered in e u , and it is swept again with a energy improvement in the inte- 
gration step AE, making a new Eq determination in which the e u and T u are recalculated. 
This recurrent relation is followed until the e u and T u values became stable. In this way we 
search these resonances, sweeping successive energy ranges (sized about 2 x 10 3 Ai? centered 
in e u ), by improving the integration energy step. The Woods-Saxon central and spin-orbit 
parameters (see TABLE HJ) are adjusted by making the stabilized energy e u reproduce the 
experimental Q value. In FIG. we present a typical example of this calculation for 147 Tm 
for which the potential was adjusted to reproduce the experimental ground-state energy. In 
TABLE [H] our results are presented in comparison with the other methods showing a good 
agreement between the data. The small differences between our calculations and the exper- 
imental results are due to the fact that we assume a pure single-particle configuration in a 
spherical model while the real system may have soft deformations or mixing configurations 
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TABLE I: The Woods-Saxon central and spin-orbit parameters used in the calculations. 
V R (r) = V 0R f(r) ; V ls (r) = V ols {^f ±/'(r)f- s 

= i +e (r-R)/ a j R = 1 - 1 ^ 1/3 ; a = 0-75/m ; V ols = 6MMeV 
Nuclide VoR(MeV) Nuclide V m (MeV) 

109 I 61.699 155 Ta 61.955 

145 Tm 63.079 W7 Ir 59.700 

147 Tm 60.398 185 5i 57.875 

151 Lu 63.000 



TABLE II: The half-lives calculated in this work are compared with the other proton emitters 
results. The Woods-Saxon and spin-orbit parameters are adjusted to fit the experimental energies. 
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Fief. [10] 
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Nuclide Orbit Q(MeV) T 1/2 


T l/2 


T l/2 


T l/2 



109j 


2c?5/2 


0.829 


7.5/j.s 


10/is 


(100 ± 5)/xs[2] 


145 Tm 


1^11/2 


1.728 


1.5/is 




(3.5 ±0.1)/xs[5j 
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2^3/2 


1.132 


280.0/xs 


210^s 206.8/zs 


(360 ±40)^2] 


151 Lu 


1^11/2 


1.255 


165.2ms 


60ms 58.4ms 


(130^°)msf2] 
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1^11/2 


1.776 


8.6/is 




(12l|W8] 


167 j r 


3^1/2 


1.086 


193.3ms 


36ms 


(110 ± 15)ms [4] 


185 Bi 


3*1/2 


1.611 


3.6/is 


3.2^s 3.17^s 


(44 ± 16)/is[3j 



altering the occupation probabilities. 

In the light proton unbound nucleus 11 N, the instability requires shorter decay times 
than those in the heaviest ones. In this nucleus the ground-state and the first excited state 
nave spin-parity inversion a8 in its mirror "B.. Recent experimental dataQ Q repot, 
the energies and widths of the ground and excited states in 11 N. In Ref. the measured 
resonances are identified by adjusting the potential to reproduce the best fit to the data. 
Using the same Woods-Saxon central and spin-orbit parameters as those given in Ref. 
we calculate the ground-state intruder | + and the other excited states. The results are 
showed in TABLE II I II in comparison with the data of Refs. Q and Q- 

We can see that 

the resonant energies present good agreement with the experimental data while the widths 
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FIG. 1: The half-life for U7 Tm is calculated by improving the integration energy step, AE, that 
means a refinement in the energy range where the resonance is evaluated. The results became 
stable when the energy precision is improved. 

have small discrepancies. The major variation in the results is due to the | + ground-state 
level which is a broad low energy resonance. In this case, the localization of the resonance 
shows up as the opposite problem to the one of the heaviest proton emitters, because here, 
the resonance width is as broad as the value of its position in energy, indicating that the 
half-life is approximately identical to its short oscillation time inside the well, which brings 
some uncertainty in the evaluation of the decay properties. 

In short, the half-lives of some proton emitters are evaluated by implementing a careful 
numerical treatment to solve Schrodinger equation in a continuum discretization approach. 
The results are in good agreement with the other theoretical and experimental data, showing 
that this method offers an efficient alternative for predicting the ground and excited states in 
proton emitting nuclei and these are necessary elements in more complex nuclear structure 
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TABLE III: The calculated single-particle resonances for unbound N are presented in comparison 
with the experimental data. 



11 N This work Ref.[l£] Ref.[14] 

r E(MeV) T(MeV) E(MeV) T(MeV) E(MeV) T(MeV) 



i + 1.38 


0.59 


1.27 


1.44 


1.63(5) 


0.4(1) 


Y 2.18 


0.51 


2.01 


0.84 


2.16(5) 


0.25(8) 


| + 3.77 


0.47 


3.75 


0.60 


3.61(5) 


0.50(8) 


| + 4.81 


0.86 


4.50 


1.27 







calculations for evaluating structures other than the single-particles ones. 
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